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Using coarse grained models we investigate the behavior of water adjacent to an ex- 
tended hydrophobic surface peppered with various fractions of hydrophihc patches of 
different sizes. We study the spatial dependence of the mean interface height, the 
solvent density fluctuations related to drying the patchy substrate, and the spatial 
dependence of interfacial fluctuations. We find that adding small uniform attractive 
interactions between the substrate and solvent cause the mean position of the interface 
to be very close to the substrate. Nevertheless, the interfacial fluctuations are large 
and spatially heterogeneous in response to the underlying patchy substrate. We discuss 
the implications of these findings to the assembly of heterogeneous surfaces. 

1 Introduction: Interfaces and Hydrophobicity 

The interaction of liquid water with oily components in aqueous solution is central 
to many phenomena pj. These phenomena are called "hydrophobic effects." The hy- 
drophobic effects leading to robust assembly, such as micelle formation and protein 
aggregation, follow from the nucleation of water- vapor- like interfaces adjacent to suffi- 
ciently extended hydrophobic surfaces [5] . This article reports on computer simulation 
studies of the behavior of such interfaces when the surfaces contain hydrophilic regions. 

Our analysis focuses on probability distributions for density fluctuations in water, 
distributions like P„ ( A^) , which stands for the probability that the centers of N water 
molecules are found in a volume v. We consider such distributions because arrange- 
ments of water molecules near hydrophobic solutes are similar to those of water near 
voids in the liquid. The connection between voids and hydrophobicity is found in Still- 
inger's [3] proposal that a water-oil interface is similar to a water-vapor interface. This 
proposal is supported by recent experiments [11 [6] . 

Hummer, Pratt and their coworkers \7] were the flrst to call attention to the behavior 
of Py{N). For small microscopic volumes in water, they used computer simulation to 
show that this probability is almost exactly Gaussian. Some theories of homogeneous 
liquids and hydrophobic effects are based upon the assumption that microscopic density 
fluctuations obey Gaussian statistics [8j , so the flnding of Ref. provides support for 
those theories. But more important is a remarkably simple yet quantitatively accurate 
theory for solvation free energies of small hydrophobic species [71 [9] . This consequence 
follows from the fact that the solvation free energy or excess chemical potential for a 
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solvent excluding volume A/iu, is given by [10] 

PAfi, = -lnP„(0), 

where 1/ (3 — k^T is Boltzmann's constant times temperature. Therefore, to the extent 
that Pv{N) is Gaussian, Afj,^ can be expressed entirely in terms of bulk water's mean 
density and mean-square density fluctuations, both of which can be determined from 
experimentally known quantities. With this solvation free energy in hand, hydration 
free energies of real apolar molecules can be estimated, using perturbation theory to 
account for solvent-solute forces beyond those of excluded volume interactions [H [TT] . 

But for large volumes v, or for small volumes in proximity to an extended hydropho- 
bic surface, the Gaussian approximation to the probability Pv{N) ceases to be accurate. 
Here, probabilities for fluctuations are not the same as those for the bulk liquid because 
the chemical potentials of liquid water and vapor differ by very little compared to k^T. 
As such, a large void in water can nucleate vapor-like configurations [3] , and vapor- like 
configurations increase the likelihood of large density fluctuations. Thus, for large v 
and N much smaller than its mean, {N)^, Pv{N) is much larger than the probability 
predicted by the Gaussian approximation. In other words, for large enough v, Pv{N) 
possesses a fat non-Gaussian tail at small values of N. Understanding the nature 
of this tail is important for estimating values of solvation free energies of sufficiently 
large clusters and extended surfaces of hydrophobic species. It is also important for 
estimating likely pathways by which water molecules are displaced during assembly of 
hydrophobic clusters. 

A theory for the fat tail or large length scale hydrophobicity must account for the 
presence of both liquid and vapor and the interface between them. The theory of Lum, 
Chandler and Weeks (LCW) [12] does so by partitioning the density field of water into 
two components. One is the field ns(r), which describes two-phase coexistence and 
is presumed to vary slowly in space (hence the subscript "s"). The other, Sp{r), is a 
Gaussian field, which is assumed to include all variation not captured by ns{r). The 
theory predicts the onset of large length scale hydrophobicity at a size of roughly 1 nm. 
Variations of ns(r) occur on length scales greater than 0.2 or 0.3 nm. The equations 
governing the behavior of Hg (r) are continuum versions of equations governing a binary 
field (i.e., the density field of a lattice gas or the equivalent spin field of an Ising model) 
on a three dimensional cubic lattice [131 HJ. To the extent that a spatially coarse 
resolution of the liquid density is of interest, a lattice gas field suffices to describe the 
fiuctuations of water. The role of Sp{r) in that case is implicit, affecting the values of 
parameters of the lattice gas Hamiltonian. We have adopted this approach in earlier 
work modeling dynamics of a hydrophobic chain in water [15j . dynamics of water in 
and around nanometer scale tubes [16j . and the dynamics of dimerization for two 
nanometer scale spheres in water |17j . Some additional justification for the approach 
can be found in theoretical analysis |18| and explicit atomistic modeling |19| . Examples 
of other applications of the lattice gas model to treat water fiuctuations near extended 
hydrophobic surfaces are found in Refs. [20| . 

The models we employ are defined in the next section. We then present results 
of this general coarse grained approach, examining density fluctuations adjacent to 
extended complex surfaces in a way that elucidates recent atomistic modeling [5T1 [55] . 
Some of our calculations focus directly on the fluid interface, i.e., a d = 2 dimensional 
manifold in d = 3 dimensional space. Our principal results demonstrate the variability 
of probabilities for dewetting fluctuations, and further show that the behaviors of these 
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probabilities are richer and more physically pertinent than might be guessed from the 
behaviors of their mean values alone. 

2 Models 

2.1 Lattice gas 

As noted, we use the lattice gas model as the discrete version of the slowly varying 
field of the LCW theory [12], with the parameterization of Ref. [18 . In particular the 
density field is taken as a binary field on a three dimensional cubic lattice with lattice 
spacing /. Positions of a lattice site are specified with r = xx + yy + zz where x, y, and 
z are Cartesian components. In units of I these components have integer values. The 
lattice site at position r has occupation number rir equal to 1 if the volume within that 
site is liquid-like and equal to if the volume within that site is vapor-like. This lattice 
gas model is capable of supporting liquid-vapor phase coexistence and its interfaces. 
The energy for a given binary density field {n^} is given by, 

Ei,{{nr}) = n,,ny,: - '^firir, (1) 

r,r' r 

where the primed summation is taken over nearest neighbor pairs of lattice sites, and 
at coexistence the chemical potential, /z, has value 3e. To match the surface tension, 
compressibility and proximity of water to coexistence, we use e = 1.51A:Br, I = 0.21 nm 
and ^ = 3e + 2.25 x 10~^fcBT. In the simulation of this model, we use cells that are 
periodically replicated in the x and y directions, while the upper z boundaries are 
populated with liquid-like lattice sites and the lower z boundaries serve to characterize 
patchy surfaces described below. 

Fluctuations across a liquid-vapor interface can be correlated over very long dis- 
tances while correlation lengths in the bulk liquid are very small. Thus, fluctuations 
of the interface can be interesting even while those in the bulk are uninteresting. For 
this reason we choose to simulate systems with large interfacial area (in the a;?/— plane) 
and relatively little bulk liquid. Taking the perspective to the extreme would disregard 
bulk entirely, as is done with an interface model, which we turn to now. 

2.2 Fluid Interface 

Here we adapt the two-dimensional model of a liquid-vapor interface described by 
Weeks [23^ . In this model, a liquid- vapor interface is treated as a two-dimensional 
manifold that is coarse grained onto a square lattice in a reference plane with lattice 
spacing approximately equal to the bulk correlation length of water, ^ = 0.42nm. This 
two-dimensional interface is periodically replicated in each of the two dimensions of the 
reference a;?/— plane. The interfacial profile is characterized by the set of height variables 
{/la}, where is the distance of the interface normal to the reference xy— plane, and 
a. — XX + yy refers to a lattice point on that plane. In this model, which we call the 
"fiuid interface model," the free interface has the energy function, 

a, a' 
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where the primed summation is taken over nearest neighbor pairs, and to be consistent 
with the surface tension of water, we use F^^ = O-I/cb^. 

2.3 Patchy substrate 

We couple the models described above to a patchy substrate. The substrate interacts 
with the solvent through excluded volume interactions as well as attractive potentials. 
The substrate lies in the z — plane, and is taken to be effectively infinite in the 
xy— plane, and it prevents the liquid (and its interface) from existing at z < 0. In 
addition to excluding volume, the substrate surface contains regions which attract the 
liquid or its interface. These are called "hydrophilic" regions. The other portions of the 
substrate, which simply exclude the liquid, are the "hydrophobic" regions. The surface 
is partitioned with the same length scale as the solvent and the state (hydrophilic or 
hydrophobic) is specified with a variable (Xa (1 or 0, respectively). For the lattice gas 
model, the interaction energy between the patchy substrate {ca} and the first (i.e., 
z = 1) layer of lattice sites is, 

AEi,{{nr}; {ffa}) = - ^(eo-a^aa + Ajn^s)- (3) 

a 

Here, na.,z refers to for r = a + zz, and a specifies x and y. The strong adhesive 
forces used to mimic the hydrophilic interactions come from the first term of Eq. [31 
The quantity A7 is the surface adhesive interaction which is included to capture small 
attractive interactions such as van der Waals attractions between oil and water. In the 
absence of hydrophilic regions (cTa = for all a), Aj{na^^i) /P is the mean field estimate 
of the difference between water-oil surface tension and water-vapor surface tension. In 
this case (all Ca set to zero), the choice A7 w 0.05e, gives an average density profile 
('T^a.z) consistent with that found experimentally for water-oil interfaces^. 

For the fluid interface model, the interaction energy between the patchy substrate 
{(Ta} and the interfacial profile is, 

ASw({/ia}; We.}) = W{{ha}) - 0(^c - ha){<y.£e + Ar^'), (4) 

a 

where Q{x) is the Heaviside step function which is equal to 1 for x > and zero 
otherwise, he is a cutoff distance equal to the bulk correlation length, the hydrophilic 
interaction energy is taken to be SksT (which amounts to about half the attractive 
energy per surface unit of e for the lattice gas), and the function W^({/ia}) exludes the 
interface from crossing the z — plane (it is infinite if any hg, < 0, and it is zero if all 
/la > 0). The surface adhesion AF is included to capture small attractive interactions 
(as above). For the case of no hydrophilic patches, the choice of AF^^ — OHk-QT 
gives an interfacial profile which is consistent with experimental observations for water 
adjacent to oil ,4J. 

2.4 Monte Carlo 

For the lattice gas model, we carry out Monte Carlo trajectories for {rir}. Acceptance 
rejection obeys detailed balance for the grand canonical ensemble with net energy 
function, 

E{{n,}, {fia}) = Ei^{{n,}) + AEi^{{n,}- {a^}). (5) 
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To specify hydrophilic patches, i.e., to specify {cth}, we tile the substrate with d x d 
squares thus creating a square lattice with lattice spacing d, where d is an integer 
multiple of /. At random, a fraction / of the d x d squares are made hydrophilic (i.e., 
Ua = 1 for each substrate lattice site in the square). The pattern so formed is then 
held fixed while the Monte Carlo trajectory is carried out for all rir's. Averages over 
these trajectories are then recorded, results of which are described in the next section. 

For the fluid interface model, we carry out Monte Carlo trajectories for {hg}. Ac- 
ceptance and rejection obeys detailed balance with the net energy function. 



Here, the hydrophilic patterns are created as above, but now with underlying lattice 
spacing ^ rather than I. Monte Carlo trajectories for h^s are performed and ana- 
lyzed with {cTa} fixed. The h^s evolve continuously, unlike the Ur's, which change 
discontinuously between and 1. 

3 Mean interfacial height for various substrates 

We have chosen to focus on surfaces which can be characterized by two parameters. 
These parameters are the overall fraction / of hydrophilic sites and the size d of the 
hydrophilic patches. The interface height of the lattice g defined such that 

/la + ^ is the smallest value of z, for z > 0, where n^^z is not zero. For example, if only 
the cell immediately adjacent to the surface at a is empty, ha = I. If none were empty, 
/la = 0. In other words, is the value of z for the occupied lattice facet closest to the 
patchy surface at a = (a;, y). Consider first the average interfacial height {h) f^d, where 
(• ■ ')f,d denotes the equilibrium average with patchy surfaces characterized by / and 
d. Specifically, 



where the first summation is over A^rop different realizations of the surface patterns, 
a refers to a specific realization, i.e., a — {ca}, consistent with / and d, and the 
second summation is over the Agurf lattice sites in the patchy surface (xy— plane). The 
quantity is the mean height of the interface over the surface site a = {x,y), for a 
specific realization. 



where the summation here is over the A^obs averaging time steps in a single Monte 
Carlo trajectory. 

In the absence of weak surface adhesive interactions, i.e., A7 — AT — 0, for some 
fraction of hydrophilic coverages, /, the mean interface height (h) f_d is non-monotonic 
in the hydrophilic patch size d. The results for these cases are shown in Fig. [T] This 
non-monotonic behavior, which is observed in both the fiuid interface model as well 
as the lattice gas model, demonstrates that for fixed surface composition / there is an 
optimal patch size d for attracting the fiuid. 

The mean interface height (h) f^d is an average over a surface which is locally hetero- 
geneous. The influence of the local surface structure on the mean (h) f^d can be visual- 
ized for a surface realization a with the spatially resolved average interface height h^"^ . 



E{{ha}, {aa}) = Ew{{ha}) + AEw{{ha}; {a^}). 



(6) 





(8) 
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Figure [2] shows /la for representative patchy surfaces with / = 0.03 and d — I, 21, and 
31. The figure shows that a hydrophilic patch affects interface fluctuations beyond the 
region immediately above the patch. In Ref. [22j . atomistic simulations were used to 
study water confined between nano-scale heterogeneous surfaces. Among their findings 
are that the solvent density within the first hydration layer over a nano-scale hydropho- 
bic region is increased significantly by the introduction of a border of hydrophilic sites. 
This effect is not symmetric with respect to exchange of hydrophilic and hydrophobic 
material as the hydration layer over a nano-scale hydrophilic region is affected very 
little by the introduction of a hydrophobic border. We have found (not shown here) 
that these results are duplicated with the lattice models used here. 

The non-monotonic behavior seen in Fig. [1] can be understood qualitatively through 
the mean interfacial profile /li"'' (shown in Fig. 2). For a surface with / — 0.03 and 
d = I in the absence of surface adhesive interactions the mean interfacial profile is 
only slightly perturbed by the underlying hydrophilic patches. Locally, regions of the 
surface with a relatively high density of hydrophilic patches lower the average height 
of the interface. There are however, many isolated hydrophilic patches for which the 
mean interface height shows little response. For d = 21, however, the hydrophilic 
patches pin the local interface, and because d is relatively small there are many such 
patches so that the distance between a patch and its neighbors is not very large. At 
d = 31 there are fewer hydrophilic patches to pin the interface and thus the hydrophobic 
domains are generally larger than when d — 21. The effect of adding a small attractive 
interaction between the surface and solute is significant as seen in Figs. [T] and [2l The 
mean interface is pulled much closer to the patchy surface when weak attractions are 
present. This result is consistent with the role of weak attractions uncovered in earlier 
work [51 1211 HH US] . The non-monotonic behavior seen in the average interfacial height 
is eliminated and the fluctuations are reduced. 

4 Density fluctuations and dewetting 

To explore the onset of drying-like phenomena we have considered a finite patchy 
substrate. The specific substrate size is 32 x 32 P for which we have studied density 
fluctuations in the adjacent lattice gas. In particular, we have computed the probability 
distribution, 'Pa(p) for various substrate realizations, a, corresponding to different 
hydrophilic fractions / and patch sizes d. The variable p refers to the density in the 
first two solvent layers above the substrate, 



where S{- ■ ■) stands for Dirac's delta function, iVg is the number of lattice sites in a 
single layer of the lattice gas (in this case Ng = 32^), and the averaging implied by 
(■ • ■)a is carried out by umbrella sampling [29j with the surface realization a fixed. This 
distribution coincides with Pv{N) discussed in the Introduction, in this case where v 
is the combined volume of the first and second layers of cells above the substrate with 
pattern a, and p times that volume is N. 




(9) 
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4.1 Without weak surface adhesive interactions, A7 = 

The inset to Fig. [3] shows Va [p) for the case of A7 = (plotted with open symbols) 
for representative surfaces with / = 0.03 and d = I, 31, and 5/. These distributions 
are very broad, showing significant fluctuations over a large range of densities. These 
broad distributions arises because the liquid- vapor interface wanders in and out of the 
observed volume. The distribution with d = I shows an additional narrow peak near 
p = 0, which corresponds to the dewetted state where the interface has completely 
pulled away from the patchy surface. 

The bimodal character of Va (p) is related to the non-monotonic behavior seen in 
Fig. [T] That is, the metastable "dry" state can emerge when the hydrophilic patch 
sizes are small, but this state is absent in the case of larger patch sizes. The loss 
of the metastability occurs when hydrophilic patch sizes are large enough to pin the 
adjacent interface. Consider the effect of excluding volume from a fluctuating liquid- 
vapor interface, as the substrate does. This imposes constraints on the configurations 
accessible to the interface, and the closer the interface is to substrate the more severe 
the constraints. Thus, for a purely hydrophobic substrate there is an entropic force that 
drives the liquid-vapor interface away from the surface, thereby drying the substrate. 
By adding hydrophilic patches, the entropic driving force is overcome by energetically 
favorable interactions between the interface and the substrate, but just barely for the 
case of / = 0.03 and d = I, as is evident in Fig. [3] 

4.2 With weak surface adhesive interactions 

The inset of Fig. [3] also shows Vaip) for the same surface parameters, / = 0.03 and 
d — I, 31, and 51, but with the interface adhesion turned on, i.e., A7 — 0.05e (plotted 
with filled symbols). With the interface pulled close to the surface the distributions 
shift to larger values of p (more liquid- like) . For d = I the addition of adhesive inter- 
actions destabilizes the "dewetted" state and as a result the distribution is unimodal 
and considerably more narrow than when A7 = 0. The main panel in Fig. [3] shows 
\og[Pa{p)] for / = 0.03, and d = ^,3^, and 51. The distributions exhibit non-Gaussian 
tails which are a characteristic of density distributions near phase coexistence [27l [28] . 
The departure from Gaussian behavior arises because density fluctuations in the pres- 
ence of a liquid-vapor interface occur through the translation of the interface. Fluc- 
tuations without the interface require compressing a nearly incompressible fluid. The 
non-Gaussian tails imply a far greater probability for large density fluctuations near 
the substrate than for those in the bulk. 

The quantity log['Pc((/9)] is proportional to the free energy to alter the solvent density 
in the volume adjacent to the substrate, and thus reflects the cost to dry the patchy 
substrate. There is a trivial dependence of the free energy difference between the p — 
and p = 1 state on the hydrophilic fraction /. For fixed /, however, the shapes of 
the distributions depend on the hydrophilic patch size d, and this feature is pertinent 
to the kinetic pathways to drying. Moving from the wet state to the dry state as a 
function of p, initially near p « 0.9, a surface with larger d incurs less of a free energy 
cost to dry than for a surface with a relatively small d. As p approaches p = 0, however, 
for the system with smaller d, fluctuations to lower p have lower free energy than for 
systems with larger d. As the solvent density, p, over the patchy substrate is reduced, 
the substrate and the hydrophilic patches dry. The solvent density over the hydrophilic 
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patches, the patches where fia = 1, is given by 



TOl = —— ^cra?la,l- (10) 

Figure|4]shows the average value of m\ for substrate reahzation a, {mi)p^a, as a function 
of p, where the average impHed by (• ■ ■)p^a is an average over configurations with fixed 
yO, and substrate reahzation a. In Fig. [4] we see that when d = I, the hydrophihc 
patches dry essentiahy uniformally as the solvent density, p, is decreased. When the 
hydrophilic patches are large, however, {mi) does not show a significant response 
to p until p is quite small, indicating that when d is large the hydrophilic patches are 
among the last regions to dry. Qualitatively, therefore, the initial stages of drying are 
easier for systems with larger d because there are relatively large hydrophobic domains 
from which to pull the interface. The final stages of drying are difficult when large 
hydrophilic patches exist, and comparatively easy for a surface with poorly pinning 
small patches. 



5 Spatial dependence of height fluctuations 

The distribution, Va [h; a) is the probability for observing interface height h over surface 
site a — {x, y) for surface pattern a. To visualize these distributions we define the free 
energy Aa(/i;a), 

^„(/j;a) = -fcBTlog , (11) 

where /li"'' is the most likely value of fta over the substrate at position a. Therefore, 
Aa{h] a) is the free energy required to move the interface at equilibrium from its most 
likely position to a height h. Figure [5] shows Aa{h]SL) for a surface a with / = 0.05 
and d — I. The weak adhesive interaction, A7 — 0.05e, causes the h — Q interfacial 
configuration to be favored by the interface. Fluctuations into h = I are thermally 
accessible to the interface {A{h; a) k^T) and for many regions, configurations with 
h — 21 (interface is 0.4 nm from the surface) are accessible through fluctuations of less 
than S/cbT. Recall that for d = 31 with weak adhesions turned on, the interface is on 
average pulled very close to the substrate (Fig. [T]), and solvent volume adjacent to the 
substrate has a more liquid-like density (Fig. [3]). But the distribution Aa{h;a) shows 
that local, fluctuations in the interface through the second layer are within the range 
of thermal fluctuations. On the other hand, the free energy to pull the interface off of 
a hydrophilic patch is quite large. These features lead to the contrasting shapes of the 
■pc<(p)'s shown in Fig. [31 as well as the response of {mi)p a to p shown in Fig. U) 



6 ImpUcations 

Many meso-scale solutes in nature are patchy. Hua et. al [30] characterized the distri- 
butions of hydrophobic and hydrophilic subunits on an assembling surface of proteins 
by coarse graining the hydrophobicity on protein surfaces over 5 x squares. The 
resulting hydrophobic surface distributions look like hydrophilic sites distributed on a 
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hydrophobic background, not unhke the surface models we are considering. The rel- 
ative fraction of hydrophilic sites depends upon the protein considered and coincide 
qualitatively to our model with / « 0.25 — 0.50 and hydrophilic patch size d « 2Z — 31. 

The protein surfaces studied in Ref. [30j and [31] are finite in size, extending over 
several nanometers. Our model surfaces are effectively infinite and thus lack a boundary 
term which may be significant in the actual kinetics of assembly. Nonetheless, for 
solutes with hydrophobic surfaces extending over about a 1 nm, the mechanism for 
assembly can involve the drying of the solvent volume between the two assembling 
proteins [2D1[3H[T7]. This drying event and subsequent solute aggregation is preempted 
by the formation of a vapor tunnel which arises when the interfaces surrounding both 
solutes come close to contact. The height fluctuations of the interface next to these 
solutes describe this phenomenon. Specifically, these interfacial fluctuations set the 
range over which solutes can broadcast their presence into the bulk solvent. Figure [6] 
shows the free energy Aa{h; a) for a substrate a with / = 0.25 and d = 31 chosen to 
resemble the assembling surfaces of the proteins in Ref. [30]. For surfaces consistent 
with these parameters, the typical range which is thermally accessible to fluctuations 
of the interface height is — 2Z (0.0 — 0.4 nm), implying that two solutes with similar 
surface distributions in water can form a vapor tunnel at separations smaller than about 
0.8 nm. Not suprisingly, therefore, Hua et. al ^50] using atomistic simuations with 
patchy hydrophobic protein pairs find spontaneous dewetting at surface separations of 
0.6nm or less in a reasonably short period of time (within 100 ps). 

The expulsion of water over buried hydrophilic patches presents a barrier to the 
water mediated assembly of patchy solutes. For an isolated patchy surface, large hy- 
drophilic patches remain "wet" until the final stages of drying. This is because drying 
these large patches is very costly, as manifest in the downward curving tails for c? = 5 
in Fig. [3] near p = 0. These downward tails do not exist for / = (not shown), which is 
consistent with the finding of Ref. [S^ that the addition of a hydrophilic site to an oth- 
erwise hydrophobic surface significantly slows drying. The free energy associated with 
fluctuations of the interface to h = 21 (0.41nm) over the hydrophilic patches are close 
to lOfcsT (Fig. [H]). With free energy barriers this large it is unlikely that the first stage 
of assembly occurs through a complete drying of the assembling surfaces. Thus, during 
assembly, water molecules over large hydrophilic regions are likely expelled in the latter 
stages of assembly through a mechanism different than those which are responsible for 
drying the hydrophobic regions. In fact, retaining some amount of buried water until 
the final stages of assembly can be advantageous. Specifically, water retained between 
assembling protein surfaces can aid in the final stages of successful aggregation [34l [S^ • 
Indeed, the grouping of hydrophilic regions can allow for much larger fluctuations of 
the interface over the substrate. This is demonstrated in Fig. [3] through the comparison 
of the distribution Va{p) for homogeneous and patchy surfaces. Along with results for 
patchy substrates. Fig. [3] has density distributions for systems with substrate layers 
of uniform attractive interactions with magnitude / (these are the "mean field" lines 
in Fig. [3]) [33] . Large fluctuations in solvent density towards the vapor- like state are 
less likely in the case of the homogenous substrate. We conclude that as geometry of 
a patchy hydrophobic surface considerably affects the mean behavior of the adjacent 
solvent, it also significantly affects the interfacial fluctuations and thus can strongly 
influences the kinetics of assembly. 
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Fig. 1 The equilibrium height, {h)f^d, of the liquid vapor interface over the sur- 
face layer as a function of patch size d for the lattice gas model (top) 
and the fluid interface model (bottom) at different hydrophilic fractions 
/. In each panel, the top three of curves are computed in the absence of 
weak water-substrate adhesive interactions, i.e., A7 = AF = 0, and the 
bottom three curves are computed with weak water-substrate adhesive 
interactions. 

Fig. 2 Mean interfacial height, ft.i"'' projected onto the a;?/— plane for hydrophilic 
fraction / = 0.03 and d = I (top left), d = 21 (top right), and d — 31 
(bottom left and right). The height h is indicated by the shading (scale 
at right)and the location of hydrophilic surface sites (substrate sites with 
fJa = 1) are shown with red circles. Panels (a), (b), and (c) correspond to 
averages in the absence of surface adhesive interactions (A7 = 0.0) and 
panel (d) corresponds to an average with surface adhesive interactions 
(A7 = 0.05e). The surfaces pictured are 64 x 64 P in size. 

Fig. 3 The probability distribution for density within the first two layers of lat- 
tice sites with weak adhesive interactions (A7 = 0.05e) over the substrate 
surface for specific realizations of patchy substrates with hydrophilic frac- 
tion / = 0.03 and patch sizes d — I, 31, and 5/. The curve label "MF" is 
the result for the mean field distribution, where the substrate is uniformly 
attracted to the solvent with an attractive strength equal to fe [33]. The 
inset compares the distributions for d = I (squares), d = 31 (circles), and 
d= hi (triangles) for substrates with (hollow symbols) and without (filled 
symbols) weak adhesive interactions. 

Fig. 4 The mean solvent density over the hydrophilic patches, (toi)p,q (see text 
for definition), as a function of the total solvent density in the first two 
solvent layers, p. Each curve is averaged over a fixed substrate realization, 
a, with / — 0.03 and weak surface adhesive interactions (A7 = 0.05e). 

Fig. 5 The spatial variation of the free energy, (3Aa {h; a) to displace the interafce 
at a from its most likely height to the indicated height h for a surface 
realization with / = 0.05 and with d = I and 31 in the presence of weak 
adhesive interactions (A7 = O.OSe). The surfaces pictured are 64 x 64 l'^ 



Fig. 6 The spatial variation of the free energy /3Aa{h; a) for h = 21, and a patch 
pattern realization with / = 0.25 and d = 31. The pattern of the patchy 
substrate {cFg} is not explicitly displayed but is evident in the pattern of 
(3Aa{h; a). The surface pictured is 64 x 64 P in size. 
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Figure 1: The equilibrium height, {h)f^d, of the liquid vapor interface over the surface 
layer as a function of patch size d for the lattice gas model (top) and the fluid interface 
model (bottom) at different hydrophilic fractions /. In each panel, the top three of 
curves are computed in the absence of weak water-substrate adhesive interactions, i.e., 
A7 = Ar = 0, and the bottom three curves are computed with weak water-substrate 
adhesive interactions. 



14 




(a) d = l,A-f = 0.00 



1 

i 


1 1 1 


1 1 


m 

i 

1 


1 

• 

1 1 ' — 


1 1 



(c) d = 31, A7 = 0.00 



T 1 1 _ I ss — I — ^-r 



J I I _ I L 



(b) d=2/,A7 = 0.00 




1.5 



0.5 



(d) d = 3Z,A7 = 0.05e 



Figure 2: Mean interfacial height, /li"'' projected onto the xy— plane for hydrophiHc 
fraction / = 0.03 and d = I (top left), d = 21 (top right), and d = 31 (bottom left 
and right). The height h is indicated by the shading (scale at right)and the location of 
hydrophilic surface sites (substrate sites with = 1) are shown with red circles. Panels 
(a), (b), and (c) correspond to averages in the absence of surface adhesive interactions 
(A7 = 0.0) and panel (d) corresponds to an average with surface adhesive interactions 
(A7 = 0.05e). The surfaces pictured are 64 x 64 l"^ in size. 
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Figure 3: The probability distribution for density within the first two layers of lattice 
sites with weak adhesive interactions (A7 = 0.05e) over the substrate surface for specific 
realizations of patchy substrates with hydrophilic fraction / = 0.03 and patch sizes 
d = I, 31, and 51. The curve label "MF" is the result for the mean field distribution, 
where the substrate is uniformly attracted to the solvent with an attractive strength 
equal to fe [33j. The inset compares the distributions for d = / (squares), d = 31 
(circles), and d — 51 (triangles) for substrates with (hollow symbols) and without 
(filled symbols) weak adhesive interactions. 
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Figure 4: The mean solvent density over the hydrophilic patches, {mi) (see text for 
definition), as a function of the total solvent density in the first two solvent layers, p. 
Each curve is averaged over a fixed substrate realization, a, with / = 0.03 and weak 
surface adhesive interactions (A7 — 0.05e). 
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{h)h = 0,d = 3l {d)h = l,d = 3l {f)h = 2l,d = 3l 



Figure 5: The spatial variation of the free energy, f3Aa{h;a.) to displace the interafce 
at a from its most likely height to the indicated height h for a surface realization 
with / = 0.05 and with d — I and 3/ in the presence of weak adhesive interactions 
(A7 = 0.05e). The surfaces pictured are 64 x 64 P in size. 
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Figure 6: The spatial variation of the free energy i3Aa{h;a) for h — 21, and a patch 
pattern realization with / = 0.25 and d = 31. The pattern of the patchy substrate 
{ca} is not explicitly displayed but is evident in the pattern oi (3Aa{h;sL). The surface 
pictured is 64 x 64 P in size. 
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